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Abstract 

We establish existence with sharp rates of decay and distance from the Chapman- 
Enskog approximation of small-amphtude quasihnear relaxation shocks in the general 
case that the profile ODE may become degenerate. Our method of analysis follows 
the general approach used by Metivier and Zumbrun in the semilinear case, based on 
Chapman-Enskog expansion and the macro-micro decomposition of Liu and Yu. In 
the quasilinear case, however, we find it necessary to apply a parameter-dependent 
Nash-Moser iteration to close the analysis, whereas, in the semilinear case, a simple 
contraction-mapping argument sufhced. 
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8 Existence for the hnearized problem [ij 
A A Nash Moser Theorem with losses [l] 

1 Introduction 

We consider the problem of existence of relaxation profiles 

(1.1) U{x,t) = U{x - St), lim U{z) = U± 

of a general relaxation system 

(1.2) Ut + A{U)U^ = Q{U), 



(1-3) U=(-), A=(^:^ M, Q- ^0 



in one spatial dimension, u € M", v € M'', where, for some smooth v^: and /, 

(1.4) q{u,v^{u)) = 0, ^a{d^q{u,v^{u))) < -0, 9 > 0, 
cr(-) denoting spectrum, and 

(1.5) (^11 Ai2) = {duf d,f). 

Here, we are thinking particularly of the case n bounded and r S> 1 arising through dis- 
cretization or moment closure approximation of the Boltzmann equation or other kinetic 
models; that is, we seek estimates and proof independent of the dimension of v. 

For fixed n, r, the existence problem has been treated in [YZilMaZT] under the additional 
assumption 

(1.6) det(^ -sl)^0 

corresponding to nondegeneracy of the traveling-wave ODE. However, as pointed out in 
[MaZ2t IMaZ3] . this assumption is unrealistic for large models, and in particular is not 
satisfied for the Boltzmann equations, for which the eigenvalues of A are constant particle 
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speeds of all values, hence cannot be uniformly satisfied for discrete velocity or moment 
closure approximations. Our goal here, therefore, is to revisit the existence problem without 
the assumption (|1.6p . 

The latter problem was treated in |MZ2| for the semilinear case, which includes discrete 
velocity approximations of Boltzmann's equations, and for Boltzmann's equation (semilinear 
but infinite-dimensional) in |MZ3| . We mention also the proof, by similar methods, of 
positivity of Boltzmann shock profiles in [LY] and the original proof, by different methods, 
of existence of Boltzmann profiles in |CN| . The new application here is to moment closure 
approximations of Boltzmann's and other kinetic equations, which are in general quasilinear. 

Our main result is to show existence with sharp rates of decay and distance from the 
Chapman-Enskog approximation of small-amplitude quasilinear relaxation shocks in the 
general case that the profile ODE may become degenerate. See Sections [2] and [3] for model 
assumptions and description of the Chapman-Enskog approximation, and Section S] for 
a statement of the main theorem. Our method of analysis, as in [MZ21 IMZ3] is based 
on Chapman-Enskog expansion and the macro- micro decomposition of [LYj . The main 
difference in this analysis from those of the previous works is that, due to a subtle loss of 
derivatives, in the quasilinear case, we find it necessary to apply Nash-Moser iteration to 
close the analysis, whereas in the semilinear case a simple contraction-mapping argument 
sufficed0 Indeed, we require a nonstandard, parameter-dependent, Nash-Moser iteration 
scheme, indexed by amplitude e — > 0, for which the linear solution operator loses not only 
derivatives but powers of e. In this, we make convenient use of a general scheme developed 
in [TZj for the treatment of such problems, which also arise in certain weakly nonlinear 
optics problems involving oscillatory solutions with large amplitudes or times of existence. 

We note that spectral stability has been shown for general small-amplitude quasilinear 
relaxation profiles in |MaZ3j . without the assumption (II. 6p . under the assumption that 
the profile exist and satisfy exponential bounds like those of the viscous case. The results 
obtained here verify that assumption, completing the analysis of [MaZ3] . Existence results 
in the absence of condition (jl.6p have been obtained in special cases in [MaZ4l IDYj by 
quite different methods (for example, center- manifold expansion near an assumed single 
degenerate point |DY] ) . However, the decay bounds as stated, though exponential, are not 
sufficiently sharp with respect to e for the needs of |MaZ3j . 

2 Model, assumptions, and the reduced system 

Taking without loss of generality s = 0, we study the traveling- wave ODE 



(2.1) 



A{U)U' 



Q{U) 





See Remark 17.41 for further discussion of this point. 
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governing solutions (jl.ip . where 

(2.3) q{u,v^{u)) = 0, ^a{d^q{u,v,{u))) < -9, 6 > 0. 
We make the standard assumption of symmetric-dissipativity [Y] : 

Assumption 2.1. (SD) There exists a smooth, symmetric and uniformly positive definite 
matrix S{U) such that 

i) for all U , S{U)A is symmetric, 

ii) for all equilibria U^f = {u,v^:{u)), dQiU^f) is nonpositive with 

(2.4) dimker 3f?5(i(5 = dimker = 

In ()2.4p and below, denotes symmetric part of the matrix M, i.e. \{M + M*). 
By the change of coordinates v ^ v — v^{u,v), we may take without loss of generality 



(2.5) v^{u,v) = 0, dQ 





dyq 



without changing either the assumed structure p.2p . (|2.ip or (since it is coordinate-independent) 
the property of symmetrizability. Note that symmetry of SdQ, together with (|2.4p . then 
implies both block-diagonal structure 

^=('o' L 

and definiteness and proper rank of 'RS22dvq- Likewise, symmetry of SA together with (12. 6p 
yields symmetry of SuAu and S'22^22 as well as 

(2.7) (511^12)^ = 522^21. 

We make the simplifying assumption ()2.5p throughout the paper. 
We make also the Kawashima assumption of genuine coupling [K]: 

Assumption 2.2. (GC) For all equilibria = (w, f*(u)), there exists no eigenvector of 
A in the kernel of dQ(U^). Equivalently, given Assumption \2.1\ (see JK^), there exists in a 
neighborhood J\f of the equilibrium manifold a skew symmetric K = K{U) such that 

(2.8) ^{KA - SdQ){U) >9>0 
for all U eM. 

Recall that the reduced, Navier-Stokes type equations obtained by Chapman-Enskog 
expansions are 

(2.9) /.(n)' = {b,{u)u'y, 
where, under the simplifying assumption ()2.5p . 

(2.10) /.(-)-/(-.0). 

b^{u)u := —Audvq ^21(^^,0). 

For the reduced system (j2.9p . symmetric-dissipativity becomes: 
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(sd) There exists s{u) symmetric positive definite such that sdf^, is symmetric and s6* 
is symmetric positive semidefinite, with dimker3?s6^= = dimkerfe^,. 

We have likewise a notion of genuine coupling [K] : 

(gc) There is no eigenvector of d/* in ker 6^, . 

We note first the following important observation of [Y] . 

Proposition 2.3 (jYj). Let (|2.ip as described above be a symmetric-dissipative system sat- 
isfying the genuine coupling condition (GC). Then, the reduced system (12. 9p is a symmetric- 
dissipative system satisfying genuine coupling condition (gc). 

Proof. Assuming without loss of generality (j2.5p . we find that s = S*!! is a symmetrizer, 
since sdf^ = ^n^n is symmetric as already observed, and s5* = —SiiAi2{S22dvq)^^ S22A21 
is definite with proper rank by the corresponding properties of S22dvQ together with (j2.7p . 
Computing that (gc) is the condition that no eigenvector of An lie in ker j42i, we see that 
(GC) and (gc) are equivalent. □ 

Besides the basic properties guaranteed by Lemma 12.31 we assume that the reduced 
system satisfy the following important additional conditions. 

Assumption 2.4. (i) The matrix 6*(u) has constant left kernel, with associated eigenpro- 
jector vr^, onto kerb,,, and (ii) The matrix := '/r*(i/^7r*(u)[kcrf)» is uniformly invertible. 

Assumption 12.41 ensures that the zero-speed profile problem for the reduced system, 

(2.11) f^u)' = {b^{u)uy, lim u{z) = u± 

or, after integration from —00 to x, 

(2.12) b4u)u' = f,{u)- f,{u±), 

may be expressed as a nondegenerate ODE in U2, coordinatizing u = (ui, M2) with ui = vr^w 
and U2 = (/ — 7r*)u |MaZ31 \Z1\ IGMWZ] . Next, we assume that the classical theory of 
weak shocks can be applied to (|2.1ip . assuming that the fiux has a genuinely nonlinear 
eigenvalue near 0: 

Assumption 2.5. In a neighborhood of a given base state uq, df^ has a simple eigenvalue 
a near zero, with a{uo) = 0, and such that the associated hyperbolic characteristic field is 
genuinely nonlinear, i.e., after a choice of orientation, Va • r{uo) < 0, where r denotes the 
eigendirection associated with a. 

Remark 2.6. Assumption [23] is standard, and is satisfied in particular for the compressible 
Navier-Stokes equations resulting from Chapman-Enskog approximation of the Boltzmann 
equation. Assumptions 12.11 and 12.21 are verified in [Yj for a wide variety of discrete kinetic 
modelsH Assumptions [27il and [231 on the reduced equations must be checked in individual 
cases. 

^ For example, both discrete kinetic models [PI| used to approximate the Boltzmann equation [PI| and 
BGK models [JXI 1N] used to approximate general hyperbolic conservation laws; see pp. 289-294 Note 
for each of these examples that the symmetrizer S is not constant, but depends nontrivially on U. 
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3 Chapman— Enskog approximation 

Integrating the first equation of (|2.ip and noting tliat f{u,v)± = f*{u±), we obtain 

(3.1) , jy 1 J \ 

A2i{u,v)u +A22{u,v)v =q{u,v). 

Taylor expanding the first equation, we obtain 

f{u,0) + Uu,0)v + 0{v^) = Mu±), 

or 

(3.2) Mu) + Mu,0)v + O{v^) = Mu±). 
Taylor expanding the second equation, we obtain 

A2i(n, 0)n' + Oi\v\\u'\) + Oi\v'\) = d^qiu,0)v + 0{\v\'^), 

or, inverting d^q, 

(3.3) v = d^q{u,0)'^A2iiu,0)u' + Oi\v\'^)+O{\v\\u'\)+Oi\v'\). 

Substituting (j3.3p into (j3.2p and rearranging, we thus obtain the approximate viscous profile 
ODE 

(3.4) h{u)u' = Mu) - /,(n±) + 0{v^) + 0(\v\\u'\) + 0{\v'\). 

Motivated by (j3.3p - (j3.4p . we define an approximate solution {uce,vce) of (|3.ip by 
choosing nc_E as a solution of 

(3.5) K{uce)uce = f*{ucE) - f*{u±), 
and vcE as the first approximation given by (13. 3p 

(3.6) VCE = C:,{uce)u'ce- 

3.0.1 Higher-order correctors 

Further expanding the second equation as 

A2i{u,{))u' + A22{u,Q)v' + 0{\v\\u'\ + \v\\v'\) = d^q{u,Q)v + 0{\v\^) 
and setting v = vce + v, u = uce, obtain 

A22{u,Qi){vce)' + 0{\v\\u'\ + \v\\v'\) = d^q{u,^)v + 0{\v\^) 
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or, inverting d^q, 

(3.7) V = d^q{u,0)-^A22{u,0)u' + 0{\vf + \v\\u'\ + \v\\v'\). 
Accordingly, we define 

(3.8) vcE,2 = vcE + d^q{ucE,Q)~^M2{ucE,Q)u'cE 

as a second-order corrector for v. Substituting vce,2 into the first equation and discarding 
the Taylor remainder as before, we obtain a second-order corrector uce,2 for u. We can 
continue this process of Chapman-Enskog expansion to all orders to obtain an approxima- 
tion 

(3.9) tj^E ■■= UCE,1 + UCE,2 + . . . , UCE,N 

to order A^, where Uce,i '■= Uce is the basic approximant at the first step. 
3.0.2 Existence and decay bounds 

Small amplitude shock profiles solutions of (13. 5p are constructed using the center manifold 
analysis of |Pe| under conditions (i)-(ii) of Assumption 12. 4t see discussion in jMaZ4] . 

Proposition 3.1. Under Assumptions \2.5\ and \2.4\ in a neighborhood of {uq,uq) in M" x 
R", there is a smooth manifold S of dimension n passing through {uo,Uo), such that for 
(u_,ti+) G S with amplitude e := — n_| > sufficiently small, and direction (u+ — u_)/e 
sufficiently close to v^uq), the zero speed shock profile equation (j3.5p has a unique (up to 
translation) solution ucE in U^. The shock profile is necessarily of Lax type.' i.e., with 
dimensions of the unstable subspace of df^{u-^) and the stable subspace of dfif{uj^) summing 
to one plus the dimension of u, that is n + 1. 

Moreover, there is 6 > and for all k there is independent of (ii_,ii+) and e, such 
that 

(3.10) \dl{ucE - n±)| < Cfce^'+^e-^^l^'l, x ^ 0. 
and, more generally, 

(3.11) \d^AncEA<Cke^^''^^e-''\^\, x ^ 0. 

We denote by 5+ the set of («_,«+) G S with amplitude e := |u+ — U- \ > sufhciently 
small and direction (ii+ — u-)/e sufficiently close to r(no) such that the profile uce exists. 
Given {u-,u+) S 5+ with associated profile uce, we define vqe by ()3.6p and 

(3.12) mE--=in^E,^^E)- 

It is an approximate solution of (j3.ip in the following sense: 
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Corollary 3.2. For fixed ti_ and amplitude e := — u-\ sufficiently small, 



K ■.= f{u^SE,v'6E)-Mu±), 

/-N -N \l /-N -N \ 

T^v ■=9{uce,vce) -q{uce,vce) 

satisfy 

where Ck,N is independent of (n_, u+) and e = \u-^- — u^\. 

Proof. For = 0, A; = 0, bounds (j3.14p follow by expansions (j3.2p and (j3.3p . definitions 
(j3.5p and (j3.6p . and bounds (j3.10p . Bounds for N, k > follow similarly. □ 



4 Statement of the main theorem 

We are now ready to state the main result. Define a base state Uq = {uo,0) and a neigh- 
borhood U = U^: X V. 

Theorem 4.1. Let Assumptions \2.I\ \2.^ and \2.4\ hold on the neighborhood U ofUo, with 
f,A,QG C°°. Then, there are Eq > and 5 > such that for u+) € 5+ with amplitude 
£ := — < ihe standing-wave equation (|2.ip has a solution U inU, with associated 
Lax-type equilibrium shock (n_,U-|-), satisfying for all k, N : 

\dl{tJ-USE)\<Cu,Ne'^''+\~'^\^\, 
(4.1) \dl{u - n±)| < Cfee'^+ie-^^l^l, x ^ 0, 

where Uce = {ucEiVce) is the approximating Chapman-Enskog profile defined in (j3.12p . 
anrf Cfc, Cyfc^jv fi^'^e independent of e. Moreover, up to translation, this solution is unique 
within a ball of radius ce about Uce in norm e~^^'^\\-\\x,2 +e~^/^ ||(9a; • ||l2 + • • ■-\-£^^^^'^ ll^^x ' IIl^, 
for c > sufficiently small and K sufficiently large. (For comparison, Uce — U± is order e 
in this norm, by H.l^i (ii)-(iii).) 

Bounds (|4.ip show that (i) the behavior of profiles is indeed well-described by the Navier- 
Stokes approximation, and (ii) profiles indeed satisfy the exponential decay rates required 
for the proof of spectral stability in |MaZ3j . From the second observation, we obtain 
immediately from the results of [MaZ3] the following stability result. 

Corollary 4.2 ( [MaZ3j ). Under the assumptions of Theorem \4-l[ the resulting profiles 
U are spectrally stable for amplitude e sufficiently small, in the sense that the linearized 
operator L := dxA(U) — dQ{U) about U has no L^ eigenvalues A with 5?A > and A 7^ 0. 
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Proof. In |MaZ3j . under the same structural conditions assumed here, it was shown that 
smah-amphtude profiles of general quasilinear relaxation systems are spectrally stable, pro- 
vided that \U'\^^ < C\U+ - |i7"(2;)| < C\U+ - |C7'(x)|, and 

(4.2) ||l + ,g„(,)ii„|<C|f/,-t/_|, - (.^.^cJ^Kw) • 

where r(no) as defined in Theorem 14.11 is the eigenvector of d/* base point Uq in the 
principal direction of the shock. These conditions are readily verified using (j4.ip . □ 

The remainder of the paper is devoted to the proof of Theorem 14. li 

5 Outline of the proof 

5.1 Linear and nonlinear perturbation equations 

Defining the perturbation variable U := U — UcE' where ^ (|3.9p . we obtain from 

(|3.1|) the nonlinear perturbation equations $^(C/) = 0, where 

(5l)^^(U)-=( h{tJhE + U)-UU-) \ 

Formally linearizing about an approximate solution U, we obtain 

(5-2) m'iu)u=( , , ^rl^r ... 

^ ' \ I \ I \A21u' + A22V' - Q22V - bU 

where 

(5.3) A = df{U'cE + U), Q22 = d,q{UljE + U), 
and 

(5.4) bU = {d{A2uA22)mE + U)U)mE + U)'- 
The associated linearized equation for a given forcing term F is 



(5.5) {^')'{U)U = F 
We have also 

(5.6) m"{U){U,U)- 

where Nj{U) are quadratic forms depending smoothly on U. 



Ni{U){U,U) 
N2{U){U,Uy + Ns{U){U,U)J ' 
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5.2 Functional analytic setting 

The coefficients and the error term TZ are smooth functions of UcE' and its derivative, so 
behave Hke smooth functions of ex. Thus, it is natural to solve the equations in spaces 
which reflect this scaling. We do not introduce explicitly the change of variables x = ex, 
but introduce norms which correspond to the usual H'^ norms in the x variable : 

(5.7) =^^11/11^, + + 

We also introduce weighted spaces and norms, which encounter for the exponential decay 
of the source and solution: introduce the notations. 

(5.8) <x>:=(x2 + 1)1/2 

For 5 > (sufficiently small), we denote by ^ the space of functions / such that e^"^^^^ f G 
equipped with the norm 

(5.9) ll/lk|,,=e^E^"'ll^''^'^^'/ll^- 

k<s 

Note that for 6 < 1, this norm is equivalent, with constants independent of e and 5, to the 
norm 

For fixed 6, introduce spaces Eg := ^ with norm || ■ \\s = \\ ■ Wh" g and := f j 



with norm I ( ) L = ||/||„s+i + 



5.3 Nash Moser iteration scheme 

Lemma 5.1. \^{0)\hi^ < j^^ allO<s<s, some C > 0. 

Proof. Immediate from p.l4p and (j5.7p . □ 

Lemma 5.2. <I>^ is Frechet differentiable from H^^^ 5, for all s > 0, e > 0, 6 > 0, 

and, for sq > 1, all s such that so + l<s + l<s, and all U,V,W H^~^^ , 

(5.10) \^'{U)\s < Co{l + \U\hs+i + \U\^s,+i\U\hiJ, 



(5.11) im'iU) . V\s < Coi\V\j,.+i + 1^1^.0+1 If/ln^+i) 

£,0 £,6 " ^ 

and 



(5.12) 



W)"{U) ■ {V,W)l < Co{\V\^y + AW\HS+l + \V\j,s+l\W\^y^ 



+ \U\j^=+l\V\„so + l\W\„so + l), 

£,S £,S 

where C is uniformly bounded for |C/|„so+i < C, for any fixed value of 6. 
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Proof. Standard, using Moser's inequality, definition (j5.2p . the fact that | ■ is a fixed 
weighted norm in coordinates x = ex, and working in x coordinates, with dx = edx- □ 



Proposition 5.3. Under the assumptions of Theorem \4-l\ there are Sq > and 6 > 
such that for all e G]0,eo], for e G]0,eo]; ^ € [0, (^o]; equation (|5.5p has a solution operator 
^^{U) (i.e., there exists a formal right inverse for {^^y{U)), such that, for all s such that 



So + 2 < s + 1 < s, So = 3, F = j G Fs, and U G H^J such that 

(5.13) \U\„so+2<Ce, 

there holds the estimate 

^'{U)F\\^,^^ < Ce-W\U\\jjs+.\Fl,+2 + \\F\\s+i) 

"^(||f^llH'S + l(||/llH- = 0+^ + Ikll H-so+a) + (||i^||r7S + 2 + IIS'II H-s + l))) 



(5.14) 



where C = C{\U\„sq+2) is a non- decreasing function of so+2. 

E,5 £,(5 

The proof of this proposition, carried out in Sections [IHS] is essentiahy identical to 
that of the corresponding proposition (Prop. 5.2) of |MZ2j in the semilinear case. Once 
it is established, existence and uniqueness follow by the abstract Nash-Moser theorems 
developed in [TZ], reproduced for completeness in Appendix lAl 

Proof of Theorem \4.1\ (Existence). The profiles Uq^ exist if e is small enough. Comparing, 
we find that Lemma [521 Proposition 15.31 and Lemma l5. II verify, respectively. Assumptions 
EH [XI and [U of Appendix [H with so = 3, 70 = 0, 7 = 1, A; = N + 2, m = r = l,r' = 0, 
and arbitrary s. Taking s sufficiently large, and applying the Nash Moser Theorem IA.4I of 
Appendix [Al we thus obtain existence of a solution of (j5.ip with |C/^|rrs+i < Ce'^^^. 

£,(5 

Defining [7^ := U^^+U^, and noting by Sobelev embedding that |/i|iTa+i controls \e^^^^^h\L^ , 

£,5 

we obtain the result. □ 

Proof of Theorem \4.1\ (Uniqueness). Applying Theorem I A. 5 1 for sq = 3, 70 = 0, 7 = 1, 
k = 3, m = r = 1, r' = 0, we obtain uniqueness in a ball of radius ce in H^q, c > 
sufficiently small, under the additional phase condition ()A.19p . We obtain unconditional 
uniqueness from this weaker version by the observation that phase condition ()A.19P may be 
achieved for any solution U = Uce + U with 

by translation in x, yielding Ua{x) := U{x + a) = Uce{x) + Ua{x) with 

Ua{x) := Uce{x + a) - Uce{x) + U{x + a) 
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so that, defining (/> := U' /\U'\, we have da{(l),Ua) ~ {(I),Uce + U') = {(p, {I + o{l))U' + U') = 
(l + o(l))|C7'| ~ and so (by the Imphcit Function Theorem apphed to h{a) := s~'^{4', Ua), 
together with the fact that {4>, Uq) = o(e) and that {(j), U'j^g) ~ \U'j^g\ ~ e^) the inner product 
(i;^, Ua), hence also Hi/a may be set to zero by appropriate choice of a = o(e~^) leaving Ua 
in the same o(e) neighborhood, by the computation Ua — Uq daU ■ a ~ o(e^^)e^. □ 

It remains to prove existence of the linearized solution operator and the linearized bounds 
(j5.14p . which tasks will be the work of the rest of the paper. We concentrate first on 
estimates, and prove the existence next, using a viscosity method. 



6 Internal and high frequency estimates 

We begin by establishing a priori estimates on solutions of the equation (j5.5p This will 
be done in two stages. In the first stage, carried out in this section, we establish energy 
estimates showing that "microscopic" , or "internal" , variables consisting of v and derivatives 
of (n, v) are controlled by and small with respect to the "macroscopic" , or "fluid" variable, 
u. In the second stage, carried out in Section [71 we estimate the macroscopic variable u by 
Chapman-Enskog approximation combined with finite-dimensional ODE techniques such 
as have been used in the study of fluid-dynamical shocks [MZH IMaZSt IPZl IZlj . 

6.1 The basic estimate 

We consider the equation 

^ ' \A21u' + A22V' + hU -Q22v) \gj 

and its differentiated form: 

(6.2) {AU' -Q + b)U = 

where h = b{UQ^)' , and A, Q, b are smooth functions of UcE + U, with ||t/|[4, [[[/(^f^H^+i 
both order e (the first by assumption, the second by estimates (j3.1ip ). We shall freely use 
below the resulting coefficient bounds 

(6.3) \d^.^^A\, [af+iQI, \d^+^K\, < Ce'^+\ \d^b\ < Ce'^+'' 

for < A; < 3 and 
(6.4) 

\di+'A\L2, \di+'Q\L^, \di+'K\L2 < Ce^+'/\e+\\U\U,), \dib\ < Ce^+^^e + \\U\Ui) 

for < j < s, where K is the Kawashima multiplier (a smooth function of ^ ). The internal 
variables are U' = {u' , v') and v. 
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Proposition 6.1. Under the assumptions of Theorem \4-l\ there are constants C, Sq > 
and 60 > such that forO < e < Eq andO < S < 60, f ^ ^eS> 9 ^ 5 O'l^dU = {u,v) ^ s 
satisfying (|6.1|) . one has 

(6-5) \\U'\\^2 +\\v\\^, <C{\\{f,f'J",g,g')h2+e\\u\\^,). 

Multiplying by symmetrizer S (block-diagonal, by assumption (j2.5p l. we obtain an ODE 
(6.6) AU' -QU + CU = F, 

where 



(6.7) A = SA, Q = SQ 

with ^Q22 negative definite, F = SF, and 




Q22 



(6.8) C = 0{u'ce)C = 0{e^)d 

comprising commutator terms and — 822^1! . 

We first prove the estimate ()6.5p for 5 = Q. Dropping hats and tildes, the ODE reads 

(6.9) AU' -QU + e^CU = F, Q = f ^ ° V 

yU W22J 

A symmetric and ^Q22 negative definite. Likewise, the genuine coupling condition still 
holds, which, by the results of [K], is equivalent to the Kawashima condition, and there is a 
smooth K = K{uce) = —K* such that '^{KA — SQ) is definite positive. Therefore, there 
is c > such that for all £ < £0 and x G M: 

(6.10) q < -eld, ^{KA - SQ) > eld. 

Lemma 6.2. There is a constant C such that for e sufficiently small, f E H^, U E H^, 
g G , and U G satisfying ()6.9p . with \\U\\2 < C£, one has 

(6.11) \\U'\\l2 + \\v\\l2 < C{\\f\\H2 + \\g\\m +4u\\l2)- 
Proof. Introduce the symmetrizer 

(6.12) S = dl + d^oK-X. 
One has 

ml o {Ad^ -Q) = ^d^oA'od^-d^oQod^- o Q' 

o K(Adcc -Q)= o ^KA od^-Re o KQ 
^{Ad, -Q) = ^A' - Q. 
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Thus 

o {Ad^ - Q) =d^ o {^AK -Q)od^ + XQ 



+ o A' o - ^AA' - oq' - o kq. 



Therefore, for U € H'^(R), (|6T0]) imphes that 

^,U\\l, + Xc\\v\\] 

-MAy\\L^{\m\\h + x\\u\\h) 



Taking 



^{SF,U)l2 > c\\d,U\\i, + Xc\\v\\i, 

- myU^lMMUh^ - \\K\\L^\\d.U\\L2\\qvh2 
-e\\C\L^\U\l. + \C'UU\f). 



A = ^||i^||i.o||(z||L= 



and using that 

(6.13) WiAYU^ + ||(Q)'!U^ = 0{e% ~ ~ .^/^(e + \\U\\,) = 0{e'/') 

and |C/|loo < = e-i/2||;7||^2 + e^/^\\U\\Hi, yields 

r'lli^ + ll^lli. < ^{SF,U)L2+e'{\\U\\l, + 
In the opposite direction, 

^{SF,U)i^2 <\\d,Uh2{\\d,{F)\\L2 + llA'llL-lli^IlL^) 

+ K\\{ny\\L4f\\L^ + \\v\\l49\\l^)- 

Using again that the derivatives of the coefficients are O(e^), this imphes that 

^{SF,U)l2 <{\\f\\H2 + \\9\\m)\\U'h2 

+ £^\\f\\L^\\u\\L2 + IbllL^bllL^, 

The estimate (|6.1ip fohows provided that e is small enough. 

This proves the lemma under the additional assumption that U € H^. When U G H^, 
the estimates follows using Friedrichs moUifiers. □ 

Proof of Proposition \6.1l This follows similarly as in the proof of Lemma 16.21 making the 
change of variables U — > e'^'^'^'lf/ and absorbing commutators. See the proof of Proposition 
6.1, pZ2] . □ 
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6.2 Higher order estimates 

Proposition 6.3. There are constants C , Cq > 0, (5o > and for all k > 2, there is C^, 
such that ^<e<eo, 6 <6o,U e H^^, U G / G and g G H^^ satisfying ^M . 

with \\U\\fj2 < Ce, there holds 

\\d',U%. +\\d'ML- ^ < C\mfJ'J",9,9')\\L-^ 
(6.14) +£''Ck{\\U'\\^k~i +e||?;||^fc-i +e||u||^2 ) 

^ ^ £,S £,S 

+ Cke'+^U\\^,+2i\\v\\Hi+e\\U\\H2). 



kTT nPlkTT — I ^xf 



Proof. Differentiating (j6.ip k times, yields 
(6.15) Ad^U - Qd^U 

wliere 

rk = -dt\{d.Q22)v)-dt\{d.A)d,U)-dt\(.d,C) U). 
Tlie estimate yields 

Wd'^u'h.^^ + iiaMi^.^^ < c{\mfJ'J",9,9')\\Li, 



for < A; < s, with = when k = 0. 

Using Moser's inequality together with (j6.3p and (|6.4p . we may estimate 

Wrkhl, < Ck{\d,Q\L^\\dt'v\\L^ + |9^QIlHI^I1l- 
+ \d,A\L^\\d',Uh2 + \d',A\MU\\L^ 

+ \d,c\L^\\dt'u\y + id'^cywh^) 

e,o e,o 

+ Cfc(e'=+^||C/||^.J|^||^.^+e'=+^||[7||^.J|C/||^.^^ 
obtaining the result by absorbing (smaller) highest-order terms from ||9^rfc|| r2 on the left- 

£,(5 

hand side. 

□ 



7 Linearized Chapman— Enskog estimate 
7.1 The approximate equations 

It remains only to estimate ||ti|t2,2^ in order to close the estimates and establish ()6.5p . To 

this end, we work with the first equation in (16. ip and estimate it by comparison with the 
Chapman-Enskog approximation (see the computations Section [3|) . 



15 



From the second equation 

A21U' + A22V' - g = dq^v, 

we find 

(7.1) V = dyq-^(^{A2i + A22di,dv^{ucE))u' + A22v'-gy 
Introducing v in the first equation, yields 

(^11 + Ai2dv^{ucE))u + A12V = /, 

thus 

{All + Ai2dv^{ucE))u' = f - Ai2v' - d'^v^{ucE){u'cE^u). 
Therefore, (|7.ip can be modified to 

(7.2) V = c^{uce)u' + r 
with 

'>' = d~^q{ucE,v^{ucE))(^A22ivy - g 

+ dv^{ucE){f' - Ai2v' - d'^v^{ucE)iucE,u))y 
This impUes that u satisfies the hnearized profile equation 

(7.3) ku' - df,u = Aur - f 
where 6* = K{uce) and df ^ := df^{ucE) = An + Audv^iucE)- 

7.2 estimates and proof of the main estimates 

Proposition 7.1. For \\U\\/i < Ce, the operator (6^,9^. — df^^){U) has a right inverse {b^dx — 
df*)^ 

(7.4) \\{hdx-df,)^h\y<Ce-'\\hh2. 

£.0 £,d 

uniquely specified by the property that the solution u = {b^dx — df*)^h satisfies 

(7.5) 4 • ^^(0) = 0. 
for certain unit vector l^. 

Proof. Standard asymptotic ODE techniques, using the gap and reduction lemmas of [MZlt 
IMaZ31 iPZl . where the assumption ||f7|[t74 < Ce gives the needed control on coefficients; 
see the proof of Proposition 7.1, [MZ2] . □ 



16 



Proposition 7.2. There are constants C, Eq > and 60 > such that for e €]0,eo], 
5 € [0,60], f G Hlg, g e H^g and U G H^g satisfying and ([73]) 

(7.6) \\U\\^, <Ce-W\f\\HS +\\g\\^, ). 

Proof. Going back now to (|7.3p . u satisfies 

kn' - df,u = 0{\v'\ + \g\ + l/'l + e>|) - /, 
If in addition u satisfies tlie condition (j7.5p tlien 

(7.7) \\u\\l2 < Ce-\\\v'\\L. + \\{f,f',g)\y +e^uh2 ). 

£ ,0 £.0 £.0 £,0 ' 

By Proposition 16.11 and Proposition 16.31 for A; = 1, we liave 

(7.8) \\U'\\^, +\Hl^ <C{\\{fJ'J",g,g')\\L2+e\\u\\^,). 

£,S £,S £,S 

(7 9) 

■ C{\\{f'j",f"',g',g")h.+e\\U'\\^, )• 

Combining these estimates, this implies 

\\v'y <C{\\if\f'\f''',g',g'')h2+e\\ifj\f'\g,g')h2+e'\\u\\^ ) 
<C{e\\{f,f',f",g,g')\\ui^^+e^u\\^,J. 

Substituting in (j7.7p . yields 

^II^^IIl^ <C(||(/,/',5)IL2 +e||(/,/',/",<?,5')llHi +^'lkllL%)- 
Hence for e small, 

(7.10) e||n||^2^^ <C(||(/,/',5)||z.^2^^+e||(/,/',r,<?,5')llHi,)- 
Plugging this estimate in (j7.8p 

(7.11) \\U'\\^, +\Hl^ +44l^ <C\\{fj',f\g,g')\\H^ 

£,S £,5 £,S 

Hence, with (17.90. one has 



\\U"h. +\\v'y < 

(7 12) 

C{\\{f',f",f"',g',g")h.+e\\{f,f',f",g,g')\\j,.). 
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Therefore, 

(7-13) \\U'\\^. +\\v\\^, +e\\u\\^, <C\\fJ',f",g,g'\\^, 

£.0 £,0 £,d £,d 

The left hand side dominates 

£,S £,S £.6 

and the right hand side is smaher than or equal to ||/||„2 + H^H f/i • The estimate ()7.6p 

£,5 £,S 

follows. □ 

Knowing a bound for ||n||r2 , Proposition 16.31 implies by induction the following final 

£,5 

result. 



Proposition 7.3. There are constants C , Sq > and 6q > and for s > 3 there is a 



constant Cg such that for e G]0,eo], ^ G [0,5o], / G 9 G -f^l^, G ^eY' 



U e satisfying ([53]), (lAlOl) . and (I5l3]) . one /las 

||f^||rrs ^ (\\U\\rTS+l\F\sf,+2 + \\F\\s + i) 

(7.14) _^ . ^'^ 

Remark 7.4. The loss of derivative on U comes from the conservative form of the lin- 
earized equations, through the microscopic energy estimates on the solution. A similar 
loss in derivative may be seen in the resolvent equation for linear hyperbolic equations in 
conservative form, XU + {A{U)uy = /; see |TZj for further discussion. We could avoid 
this by writing the differentiated equations in quasilinear form, but this would prevent us 
from integrating back to carry out linearized Chapman-Enskog estimates. That is, the loss 
of derivatives is due to a subtle incompatibility between the integrated form needed for 
linearized Chapman-Enskog estimates and the nonconservative (quasilinear) form needed 
for optimal energy estimates with no loss of derivative. 



8 Existence for the linearized problem 

To complete the proof of Proposition 15.31 it remains to demonstrate existence for the lin- 
earized problem. This can be carried out as in |MZ2j by the vanishing viscosity method, with 
viscosity coefficient rj > 0, obtaining existence for each positive r] by standard boundary- 
value theory, and noting that our previous A Priori bounds ()7.14p persist under regulariza- 
tion for sufficiently small viscosity r/ > 0, so that we can obtain a weak solution in the limit 
by extracting a weakly convergent subsequence. We omit these details, referring the reader 
to Section 8, [MZ2| . The asserted estimates then follow in the limit by continuity. 
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A A Nash— Moser Theorem with losses 

For completeness, we give in this appendix the parameter-dependent Nash-Moser theory 
developed in [TZj . specialized for clarity to the present, Hilbert space, setting. The main 
novelty of this treatment is to allow losses of powers of the parameter e — > in the linearized 
solution operator. For a proof of this result, see |TZj : for a more general discussion of Nash- 
Moser iteration methods, see [Hi \AG\ IXSR| . and references therein. 

Consider two families of Banach spaces {-E^lsesR, {-^^sIsgk, and a family of equations 

(A.l) = 0, u'eEs, 

indexed by e E (0, 1), where for all e, 

(A.2) G C^Es, Fs-m), for all s < s, 

for some m > and some s G 3ft. 

Let \ ■ \s denote the norm in Eg and || • ||<j denote the norm in Fs- The norms \ ■ \s and 
II • \\s may be e-dependent (as in our application here). We assume that the embeddings 

(A.3) Es'^Es, Fs'^Fs, s<s', 

hold, and have norms less than one: 

(A.4) I • Is < I • Is', II • lis < II • lis', s < s'. 

We assume the interpolation propertjH: 

(A.5) I • Is+a < I • Is^l • L^^M 0<a<a'. 

We assume in addition the existence of a family of regularizing operators 

Se '■ Es ^ Eg, 9 > 0, 

such that for all s < s,' , 

(A.6) \Seu-u\s<e''''\u\s>. 

(A.7) \Seu\s' <e''-'\u\s. 

Assumption A.l. For some sq G 5R, some 70 > 0, for all s such that 

SQ + m<s + m<s, 



^In (|A.5|) and below, \u\s < \v\g, stands for |itjs < C|w|s/, for some C > depending on s and s' but not 
on e, nor on u and v. 
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for all u,v,w G Eg+m, 

(A.8) \\<^>'iu)\\,<Co{l + \u\s+m + \u\ 

(A.9) ||(^>^)'(n) • v\\s < Coi\v\s+m + \v\so+m\u\s+m.), 

and 

/ . -I^s ll(^')"(^) • (.V,w)\\s < Co{\v\so+m\w\s+m + \v\s+m\w\so+m 

^ II II II \ 

where Cq = Co(e, satisfies 

(A. 11) sup sup Co < +CXD. 

Assumption A. 2. For some 7 > 0, r > 0, r' > 0, /or a// s such that 
(A. 12) sq + m + max(r, r') < s + max(r, r') < s, 

for all u € -E^+r such that 

(A.13) klso+m < e'^, 

i/ie map ($^)'(u) : i?s+m ^ -^s /^cis a right inverse ^'^(u) : 

($^)'(u)^'^(7z) =Id: Fs^Fs, 

satisfying, for all (f) G F^^^'^ 

(A.14) |^^(n)0|. < €-1^(11011 

where C is a non- decreasing function of its arguments s and \u\sQ-\-m+r- 
Assumption A. 3. There holds the bound 

(A.15) ll^^(0)||s<e', 
for some k and s satisfying 

(A.16) max(2, 1 + 70,1 + 7) < fc, 

(A.17) Cik) <s - so-m, 

where C{k) is a certain positive function (see \TZ^ ) and s & [sq + m,s — C{k)]. 
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Theorem A. 4 (Existence). Under Assumptions ETTl IA.2I and IA.3[ for e small enough, 
there exists a real sequence 9^, satisfying 6^ — > +oo as j — > +oo and e is held fixed, such 
that the sequence 



Un 



0, n^.+i := + Sg.v], v] := -^'^(7xp$^(np, 



is well defined and converges, as j ^ oo and e is held fixed, to a solution of (|A.ip in 
s + m norm, which satisfies the bound 

(A18) \u%<e''-\ 

Theorem A. 5 (Uniqueness). Under Assumptions \AA\ IA.2I and IA.31 for e small enough. 



if {^^y is invertible, i.e., is also a left inverse, then the solution described in Thm A. 4 
IS unique m a ball of radius o(e"^^^(i''ro,7)) in SQ + 2m + r' norm. More generally, if vf" is 
a second solution within this ball, then {if — n*^) is approximately tangent to Ker(<I>'^)'(ti^), 
in the sense that its distance in sq norm from Ker($^)'(u^) is o{\u'' — u^Isq)- -^'^ particular, 
z/ Ker(<I>'^)'(ii^) is finite- dimensional, then u is the unique solution in the ball satisfying the 
additional "phase condition" 

(A.19) nKer($.y(„.)(^'-^') = 0, 

where ^Y>.eT{'^'')' [w^) '^^^^ uniformly bounded projection onto Ker($'^)'(n'^) (in a Hilbert 
space, any orthogonal projection onto Ker($^)'(u^)j. 
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